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ABSTRACT 

Diffusion measurements were made at the Round Hill Field Station 
Nebraska in the year 1957. These measurements were used %G calculate 
04 , the LaGrangian standard deviation in the direction lateral to the 
mean wind for a series of experimental runs under varying conditions of 
thermal stability, and horizontal and vertical wind shear. A relationship 
between Ty and the downwind distance was derived from Sutton’s equation 
for an elevated point source. Tay was calculated for the various runs of 
data and correlated with the crosswind and downwind distances. Correla- 
tions were worked out between (i) Oy and the standard deviation of the 
wind angle {ii) Vay and the contour profile number p, and (iii) 14 and 
the stability ratio. A marked correlation was found for Tn the azimuth 
wind variation. 

The writer wishes to express her appreciation for the assistance 
and encouragement given her by Professor Frank L. Martin of the U. S. 


Naval Postgraduate School in this investigation. 
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SYMBOLS 
Generalized diffusion coefficients 
Concentration of aerosol along main axis 
Concentration not along main axis 
Height of the source above the ground 
Source strength 
Correlation coefficient 


Lagrangian autocorrelation coefficient of velocity 
component of lag 


Wind in the x-direction 

Standard deviation of azimuth wind angle 
Lagrangian deviation along y-axis 
Lagrangian deviation along z-axis 


Wind speed at 2 meters 





‘ae Introduction 

During the fall of 1957, a series of diffusion experiments from a 
steady state point source was conducted at the Round Hill Field Station 
in Nebraska. The sampling array consisted of three overlapping, inde- 
pendently operated networks at travel distances of 50, 100, and 200 
meters. During the experiments, time-mean concentrations for sampling 
intervals of .5, 3, and 10 minutes were obtained for each travel distance. 
A schematic diagram of the field installation is shown by Figure 1. The 
ten-minute network comprised of individual stations located at a height 
of 1.5 meters and were spaced at 3degree intervals along 180 degrees of 
arc: sampling stations for the three-minute network were at a height of 
1.5 meters and spaced at 1.5 degree intervals along an arc of 150 degrees. 

A sulfur dioxide generator was used to supply the tracer. A point 
source of strength 100 g per sec was required during conditions of thermal 
instability while an emission rate of half that amount was sufficient under 
nighttime thermal conditions of thermal stability. Prior to the start of the 
experiment, the tracer was permitted to traverse the entire network. The 
three sampling networks were then turned on simultaneously and each 
operated for the appropriate length of time. Aspiration of the impingers 
was provided by ten vacuum tanks. Error introduced into the concentration 
measurements by the operation of the vacuum system was estimated to be 
less than five percent. 

Meteorological instrumentation included: a cup anemometer and 


sensitive azimuth vane located at a height of two meters near the source; 
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cup anemometers and ventilated thermocouples at heights of 1.5, 3, 6, 
and 12 meters on a portable tower; and for most experiments, five bivanes 
equipped with heated thermocouple anemometers. The operation of all 
meteorological instrumentation was controlled by a timer located within 
the recording truck. A twenty-minute observation period centered on the 
ten-minute gas sampling period was employed for the meteorological! 
instruments. 
Zi Theoretical considerations: 

Assuming that the mean concentrations within a diffusion cloud 
are distributed laterally by a two-dimensional normal distribution, Sutton 
has shown that the steady-state concentration from a point source is ae 
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for a ground source. For an elevated source of height, h, the correspond- 
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Sutton‘'s reasoning, in addition to the assumption of the Gaussian distri- 
bution of pollutant, is that the concentration shall satisfy the following 
boundary conditions: 


(de C—>0, 0% bus +0 


(2). The downward transport at the earth's surface is zero. 
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(3). The rate of the total transport of contaminant through any 

downwind plane is constant and equal to the generation rate at 

the source. 

In this study expression for y and Wes of the form 

G- 
Ty = Bx ; 4p >? O 

will be sought. 

Sutton EY made certain assumptions regarding the form of the 
Lagrangian autocorrelation function R (% which relates the y-component 
of velocity with that existing for the same particle at time %y later. He 


then obtained: 
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Here Cy Co are called the generalized diffusion coefficients. 
Sutton a / has obtained formulas Cy: Co in terms of the gustiness of 
the wind. The particular formulas for Cy. C, may be found on page 251 
er /3 /. 

The expressions for (i 73 of (4), (5) are usually introduced 
into (1) or (2) at this point. Barad and Haugen jl 1 /have tested Sutton's 
theory using the equation modified in this way. However, the values they 
found did not agree with n = 2p/p+ 1. Thus in this paper, the more general 
Sutton theory of Equations (1) and (2) has been used, with the object of 


determining empirical information regarding Ty and i 
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Let oe be the axial concentration at a given downwind distance x, 
and Cy the concentration at (x,y,z) where, in the major part of the Round 
Hill data, z= 1.5 meters. In order to obtain Vy G3 as a function of x, 
the assumption was made that these parameters did not vary over the 
measuring range lateral to the axis. Therefore, using Co with y = 0 in 


Equation (2) we obtain: 


@) St = exp (- 4y* /2.04*) 


On taking logs of both sides of Equation (6) 


 Ty'= Zink*f,) 


This was the main working equation of this paper. A similar procedure 
will be outlined for (Vo , although no detailed computations for Jz. have 
been obtained. 

S Research objectives and results: 

The first problem of this study was to determine a relationship 
between Vy and the downwind distance x along the mean wind. An 
equation of the form b 
T= Bx’, bo 
was tested. Secondly, it was desired to determine the empirical relation- 
ship between Vy and x and y, where y is the distance lateral to the wind 
axis. 

A chemical cloud may be pictured ideally as shown in Figure 2. 

The data Van presented concentrations at only three downwind distances 
along the axis: at 50 meters, 100 meters and 200 meters. The concentration, 
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Figure 2. Sulfur dioxide cloud from point source 





Cg, at other points along the x-axis is also needed. For each of the 
10 runs, the known concentrations were plotted on log-log paper and c, 
for each value of x in the range 50 to 200 meters was determined. 

Under extremely stable conditions the cloud, as might have been 
predicted, showed limited spreading and very little could be deduced as 
to the relationship of Ty tox andy. In unstable conditions, the cloud 
spread rapidly vertically and laterally. Results of the calculations are 
shown in the Appendix. 

In Sutton's hypothesis, a cloud with its maximum concentration 
along the mean wind axis was assumed to be symmetrical. It was found, 
in practice, that the cloud was asymmetrical, and that maximum concen- 
tration did not quite coincide with the direction of the maximum wind 
frequency. The x-axis was therefore assumed to be along the line of 
maximum concentration. There were also other directions along which 
the concentrations reached secondary maxima, with rather closely asso- 
ciated secondary maxima of wind frequency direction. These distributions 
represented superimposed bimodal or trimodal Gaussian distributions. The 
variance V4 of the Sutton theory was assumed, however,to be that of a 
unimodal Gaussian distribution, Eq. (2), but actually varied markedly 
from one side of the major axis of the cloud to the other, and increased 
also with the length of the averaging time. 

There was,in most instances, a rapid decrease of 14 near the axis. 


For a possible explanation, let us examine the equation: 
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For points near the axis, thatis, for c, very nearly equal to Co: the ratio 
C)/c, approaches one. The natural log of one is zero, and values of the 
log function near one are very small. The concentration C) near the axis 
does not decrease rapidly at first as y increases; however y increases to 
the second power and the value of 4 increases very rapidly. Ata short 
distance from the axis the concentration Cy begins to decrease appreciably 
and the value of Jy then begins to decrease with an increase in y. This 
peaking of ly also happens to some extent aS secondary maxima concen- 
tration-directions are sampled. 

Noting that there was a straight-line decrease ny with a decrease 
of x, when these coordinates were plotted on log-log paper, for impingers 
well-removed from concentration maxima, an empirical relationship be- 


tween a and xbym was sought. The following relationship was tested: 
fw 
7) T4 - Dy “4 


where D, m and b are unknowns. The results from use of Eq. (7) indicated 
that both m and b were negative. However, observationally b should be 
positive when x is varied, while keeping y constant. Hence the assumed 
dependence upon y was rejected. 


The equation was modified to the following form: 


(8) Tay = Bx* 


In order to get a good fit for b in each run, the statistical method of 
least Squares was used. 


Eq. (8) was reduced to the following form: 
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Summing up this equation into the least squares form: 
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log x 


(10) 2 log Te ae loqB B logy + bE (lox), 


where N is the number of points for which corresponding Ty and x are 


available on the measuring arc. 


From these equations, b and B were 


obtained from each ten-minute run except #1 where the thermal stability 


» prevented the cloud from spreading to any extent. 


sented together with those of OA 


meters will be defined on page 10 in Table l. 
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Values of B and b and micrometeorological parameters. 


The values of Uy which appear in Table ] wereobtained from Vor, » The 


stability ratio, SR, was defined from the equation: 
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which is essentially a measure of the Richardson number. The profile 
contour number p, which was also used extensively in the Sutton theory, 


is defined by the equation: 


@)* wow a 


Values of winds at 12 and 6 meters from reference A were used in 
computing p. 
4. The correlation analysis. 

The parameters B and b which specify (ag were linearly correlated 
with each of t.. SR and p. The number of independent cases in each 


sample was N=9. The results of the correlation were as follows: 


1 rB y) = -.457 
2 r(b,—,) = -.844 
2 (Bye) — — 207, 
4 r(b,SR) = -.0946 
5 r(B,p) = B= 47 
6 - r(b,p) = -.601 


The celpheatirn correlation coefficient which is significant at the 
.05 level, based on nine independent cases is .816; at the .10 level, 
the minimum is .796. 

Correlation (2) is significant at the .05 level. 

Correlations (1), (5) and (6) may be significant but not at the 
levels of .05 or .1. The ctners were not significant at all. 
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It can be concluded that it may be possible to specify ly = Bee on 
the basis of the micrometeorological parameters , and p although a larger 
sample of diffusion data would be necessary to establish adequate signi- 
ficance levels. 

Sr Suggestions for future research. 

One of the assumptions made was thatTy did not vary appreciably 
over a lateral cross-section (at fixed x) through the cloud. The value 
obtained using this assumption was therefore a first approximation. It 
was found that the values thus obtained did vary laterally especially on 
passing through a secondary maximum of concentration. It would be 
desirable to determine values of Vy which are relatively constant with 
y in spite of secondary maxima of wind frequencies. 

The data for this experiment can also be analyzed for Vz : 

Making the assumption, again, that Ty and do not vary appreciably in 
the measuring range lateral to the AX1S, the following equation was derived 
in the same manner as the equation for , uSing the levels of 1.5 and 


2.9 meters for which data are ( (ho) “ ( by 6 ke 
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From this equation, A may be found, and from that Vz 

In order to establish firmer tests of correlation significance, a 
greater number of runs should be analyzed using Eq. (8). 

It is also desirable to test the diffusion model of Cramer et al LUE 
in order to see if their diffusion parameters are more nearly predictable 


from the meteorological variables than of this study. 
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Right side - 10 minutes 


Run 6 — 200 meters (continued) 


x y e4 Co TH 
190.2 61.6 8.43 29.3 32 
186.7 7.7 0.662 30.0 25.8 
Sz .7 31.3 0.200 31.0 20 oO 


Run 7 - 100 meters 


99.86 &.23 35.8 3069 13.8 
99.45 10.45 35.3 38.0 23 
wo.77 15.64 34.8 39.9 29.4 
oe 82 20.29 26.1 40.0 Zand 
96.59 25.88 15.4 41.0 18.5 
oe,,11 30.90 8.05 41.2 16.8 
93.36 35.84 4.66 41.5 17.1 
91.36 40.67 2.21 41.8 16.8 
89.10 45.40 0.500 Az 3 15.2 
56.60 50.00 0.065 43.0 13.8 


Run 7 - 200 meters 


189.7 10.5 12 0 i.” 20.6 
re . 9 2029 8.38 13.8 20:9 
197.5 31.6 5.41 14.0 22 
195.6 41.6 3.00 eZ 23.6 
He 2 51.8 1.33 14.5 Fags | 
190.2 Gis 0.615 mS .0 24.4 
186.7 (Ales 0.239 bere 24.8 
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Right side - 10 minutes 


Run 8 - 100 meters 
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Run 8 - 200 meters 
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Run 9 - 100 meters 
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Right side - 10 minutes 
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Run 9 - 100 meters (continued) 
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Run 9 - 200 meters 
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31.3 


41.6 
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Run 19 - 100 meters 
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Right side - 10 minutes 


Run 10 - 100 meters (continued) 


y Cy - 
54.46 14.2 4S. 
S828 10.3 a. 
62.93 4.87 49. 
66.91 29 SUF 
“On7 | 1.94 SSH 
72.31 ‘hele GH. 
7 2 An 2 aS 
80.90 3.50 ILE 


Run 10 - 200 meters 


10.5 3.03 4 
ZA, 9 3.79 4 
31.3 4.26 4. 

4.6 3.23 a. 
ob. 8 Lieae che 
61.6 1.46 D. 
7 1.00 oe 
oo. 3 1.43 Oe 
90.8 yar (iS) as 
100.0 3.49 6. 
Me. 9 2.41 6. 
117.6 ol (e 
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Right side - 10 minutes 


Run 10 - 200 meters (continued) 


x x Cl Co 4 
Poa o4 LZong 2.78 oie) 85.0 
148 .6 133.8 A Ole: SS. 63.8 
141.4 141.4 0.597 10.0 yell) 
133.8 148 .6 ORS 2/ 11.0 60.3 
Zor. Poo 34 0.967 13.0 68 .6 
7; 6 161.8 Le? 17.0 O97, 
108 .9 O77 0.166 19.0 Sad 
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